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OS. 

Q\ , We study the effect of a background flux string on the vacuum energy of massive Dirac fermions 

in 2+1 dimensions confined to a finite spatial region through MIT boundary conditions. We treat 
two admissible self-adjoint extensions of the Hamiltonian and compare the results. In particular, 
for one of these extensions, the Casimir energy turns out to be discontinuous at integer values of 
the flux. 
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I. INTRODUCTION 



In a by now classic paper, Y. Aharonov and D. Bohm Q| pointed out that a locally trivial vector potential can, 
however, give rise to observable effects in a nontrivial topology. Since then, the relevance of Aharonov-Bohm scenarios 
both in particle physics and condensed matter was recognized. 

More recently, much attention has been paid to the inclusion of spin, mainly in connection with the interaction 
of cosmic strings with matter [2H§. In this context, the need to consider self-adjoint extensions of the radial Dirac 
Hamiltonian was understood |7 1| . 

Thereafter, the vacuum properties of Dirac fields in the background of a singular p-[L!f| as well as extended [|15|Jl6[ 
magnetic flux string have been extensively studied. Scalar field and the electromagnetic field have been considered as 
well; some pertinent references are fl7|-|20|] . 

As is well known, the presence of background fields modifies the energy spectrum, thus giving rise to a nontrivial 
vacuum, or Casimir, energy fl2ll,E2j. Furthermore, the Casimir energy is altered by the presence of boundaries, and 
the consequent imposition of boundary conditions on the quantum fields. For Dirac fields, many examples of both 



5_i ■ situations have been studied in the literature (see, for instance, [|16],[23 25 1). 

In particular, the combined effect of a classical magnetic fluxon and MIT boundary conditions on the vacuum energy 
of a massless Dirac field in 2+1 dimensions was treated in Efl]. There, one of the possible self-adjoint extensions of 
the radial Hamiltonian was considered. 

In this paper, we will consider the more realistic case of massive fermions. In this context, it is important to 
mention that, despite general belief, in the presence of curved boundaries the effect of a mass is not exponentially 
small p7| as it is for parallel plates [E2[. On the contrary, in some situations it might even lead to a sign change in 
the Casimir force and is by no means negligible P7[ . In the present context, we will see that properties like existence 
of a minimum or continuity of the Casimir energy as a function of the flux depend on the mass which underlines 
its crucial importance. Moreover, we will analyze two possible self-adjoint extensions, both known to be compatible 
with the presence of a Dirac delta magnetic field at the origin [£8| . The resulting Casimir energies are quite different 
showing that the one-parameter family of self-adjoint extensions describes nontrivial physics in the core, see also 0. 
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The organization of the paper is as follows. 

In Sec. Owe summarize the generalities of the model. 



In Sec. Ill, we present a discussion of self-adjoint extensions distinguished by the behaviour of the wave function at 
the origin. We determine the energy cigenfunctions corresponding to two different cases of these extensions. The first 
one is a minimal divergence extension. As shown in p3,Q4|, it arises when imposing Atiyah-Patodi-Singer boundary 
conditions p9|-|34| at a finite radius, then shrinked to zero (an idea first suggested in |3£j]). The second one follows 
from the zero radius limit of a cylindrical flux shell [Q,|]{]] ■ 



In Sec. IV, the implicit equations for the energy spectrum are found in both cases, once the theory is confined to 
a circle of radius R and MIT conditions (see, for instance, p6[ and references therein) are imposed at the exterior 
boundary. The expression of the Casimir energy for both types of behaviour at the origin is given in the framework 
of the £- function regularization [pTJ-Eof . 

In Sec. [V], we evaluate the vacuum energies, following the methods developed in [ p6pl| - |4^p7[ | . 



Finally, Sec. VI contains a discussion of the results. 



II. SETTING OF THE PROBLEM 

We study the Dirac equation for a massive particle in 2 + 1 dimensional Minkowski space 

(i ft-4-m)^ = (1) 

in the presence of a flux string located at the origin, i.e., 

H = V A A = -S(r)e z , (2) 

r 

where k = ■$- is the reduced flux. 

It: 

We assume the flux string to be radially symmetric; so, a gauge can be chosen such that the (covariant) vector 
potential is given by 

Ae(r)=--, for r>0. (3) 

r 

We will consider the chiral representation for the Dirac matrices 

7 =P3®cr3, 7 1 =ip 3 ®(72, l* = ~IPZ ® &1 , (4) 



which, together with 



7 = ip2 <8> (T , (5) 



give a closed Clifford algebra. 

Then, the eigenvalue equation for the Dirac Hamiltonian takes the form 



where the two-by-two blocks are given by 



"J h_ | »£=£»,_■. »i) 



»=it£). 



and we have introduced L = -ie ie {-d r + B), L+ = ie- w (d r + B), B = -^d e - f. 

(Notice that these two "polarizations" correspond to the two inequivalcnt two by two irreducible representations of 
the gamma matrices |3|). 

The general solution to Eq. (m can be written as a combination of 

< = (*f). •s°-(; ; ). m 

with 



H±1>% 



E^>% 



(9) 



After separating variables, and for nonintcger n = £ + a (where i is the integer part of the reduced flux, and a its 
fractionary part), the eigenfunctions in Eq. (W) turn out to be 



^e (r, l 



{EZ-ooti(r)e inf> \ 



( ££=_«, {K J n -K (kr) + B+ J K _ n (kr)) e* 



EOO 
n— - 



At (A+ J„ + i_« (fcr) - £+ J K _„_! (kr)) e 



i(n+l)6 



E"=-oo (^n Jn- K (kr) + B- J K _„ (At)) e* 



v r°° 

\ Z_- /n— — c 



(10) 



(A- J n+1 _ K (fcr) - B- J«_„_! (fcr)) e l (" +1 ) e / 



where fc = +-^/(£' 2 — to 2 ). 

(Of course, for integer n, a linear combination of Bessel and Neumann functions must be taken). 



III. BEHAVIOUR AT THE ORIGIN 



As is well known [HJ3J3I , the radial Dirac Hamiltonian in the background of an Aharonov-Bohm gauge field requires 
a self-adjoint extension for the critical subspace n — L In fact, imposing regularity of all components of the Dirac field 
at the origin is too strong a requirement, except for integer flux. Rather, one has to apply the theory of Von Neumann 
deficiency indices J44JJ , which leads to a one-parameter family of allowed boundary conditions |q] , characterized by 



i lim (mr) a Qe(r) sin 1 I = lim (mr) a fe(r) cos 1 

r^o v ' ye v ' \ 4 2 r^a v ' Jl K ' \ 4 2 



(11) 



Here, ± parameterize the admissible self-adjoint extensions of H± respectively. Which of these extensions to choose 
depends on the physical situation under study. 

Throughout this paper we will, for non-integer n, consider two different behaviours at the origin. The first one, 
from now on called behaviour I, is characterized by 



e* = 



for a > | 
for a < % 



(12) 



As shown in Rcfs. 23 24f|. this is the extension arising when boundary conditions of the Atiyah-Patodi-Singer (APS) 
type ]2^-|3^| are imposed at a finite radius, then taken to zero. 

The second self-adjoint extension we will consider, from now on called behaviour II, corresponds to 



e ; 



-f for K >0 
f for K < 



(13) 



As shown in Ref . jq] , this extension arises when a finite radius flux tube is considered, thus asking for continuity of 
the components of the spinor, and then shrinking the radius to zero. 

Outside the critical subspace, the eigenfunctions in Eq. ( |l0| ) are determined by the requirement of square integra- 
bility at the origin, and are thus identical for the behaviours I and II. They are given by 



B+J, 



(K 



£-\-a — n 



(kr) e l 



n-(r,o) 



,i(n+l)8 



*E^ B n J l+a-n-l(kr)e 

B-J e+a _ n (kr)e m0 

V*^^n^+a-™-l(fcr)^ (Tl+1)e 



\ 



(14) 



and 



* 



n>i+l 

E 



(r,9) 



( A+J n -i-a(kr)e me \ 

A-J n _ t _ a {k,r)e m9 



\-*E^A-J n+1 ^ a {kr)e^) s 



(15) 



In the critical subspace (n = £), the eigenfunction for behaviour I is given by 



¥ E {r 



i{i+i)e 



( B+J a {kr)e m 
B^J a (kr)e m 



\ 



for a > 



(16) 



and 



^ E (r,t 



(A+J^ a {kr)e m \ 

-i^A+J^ a (kr)e^+^ 
AjJ- a (kr)e m 



for a < 



(17) 



It is easy to see that this extension satisfies the condition of minimal irregularity (the radial functions diverge as 
r — ► at most as r~~ p , with £><§)• Moreover, it is compatible with periodicity in k, a natural requirement when the 
origin is an excluded point. 

When behaviour II is imposed at the origin, the eigenfunctions in the critical subspace are given by Eq. (|l6|) for 
K > 0, and by Eq. (jl7|) for k < 0. 

It is worth pointing out that, for integer k = £, both APS boundary conditions and the finite radius flux tube lead, 
when taking the singular limit, to the requirement of regularity of all components at the origin. In this case 



/ A n Jn—K, (Ax J e 



®E{r, 



E 



n— — oo 



' E-m 

A n Jn—K ykr) e 



A+J n+1 _ K (kr)e« n +W 



\ 



(18) 



\^^A-J n+1 ^(kr)e^+^ J 



IV. THE THEORY IN A BOUNDED REGION. ENERGY SPECTRUM AND CASIMIR ENERGY 



From now on, we will confine the Dirac fields inside a bounded region, by introducing a boundary at r = R, and 
imposing there MIT bag boundary conditions. 
The Casimir energy is formally given by 



Er 



\E>m E<-m / 



(19) 



where p represents all indices appearing in the eigenvalue equation that arises after local MIT conditions are imposed. 
In doing so, one must consider a boundary operator B which, with the representation of the Dirac matrices given in 
Eq. (0), is also block-diagonal, and can be written as 



B = l-i 

where n is the exterior normal, and 



1 + j(7 i n i +7 2 n , ) 



B+ 02x2 
02x2 B— 



(20) 



B± 



1 ±ie~ l ° 
Tie 10 1 



(21) 



Consider, in the first place, behaviour I (II) for a > 1/2 (k > 0). Then, the eigenvalue equations for the upper (+) 
polarization are 



k 

J n+a (kR) = — J n -i +a (kR) for n = 1, ..., oo, (22) 

h, — m 

— k 

J n -a(kR) = -= J n+ i- a (kR) for n = 1, ..., oo , (23) 

£/ — m 



coming from non-critical subspaces, and 



J a (kR) = — — Ja-i(kR) , (24) 

A — m 



which comes from the critical one. 

The eigenvalue equations corresponding to the lower (— ) polarization are 



— k 

J n+a {kR) = -— J n -i+a(kR) for n — 1, ..., oo, (25) 

hi + m 

k 

Jn-a(kR) = -=-, Jn+l-a(kR) for U = 1, ..., OO , (26) 

hi + m 



from non-critical subspaces, and 



J a (kR) = — — Ja-i(kR) , (27) 

hi + m 



which comes from the critical one. 

For a < 1/2 (k < 0), the contributions from non-critical subspaces are the same, while those due to n = £ are 



for the upper polarization, and 



J- a (kR) = —Jt-J^kR) , (28) 

h — m 



J-a(kR) - tti— Jl-„(fcU) , (29) 

h, + m 



for the lower one. 

It is easy to verify that positive energies coming from one polarization correspond to negative energies coming from 
the other. Thus, both polarizations give identical contributions to the Casimir energy in Eq. (p9|). 



Then, the formal expression for the Casimir energy (19) is 



Ec = ~2j2(k 2 + m *) 1/2 , (30) 

where the fc's are solutions of 



2 

k 



Jn+a(kR)-Jn-l + a(kR)-^Jn+a(kR)Jn-l + a(kR)=0 for U = 0, ..., OO, (31a) 

Jl_ a {kR)-Jl +1 _ a {kR) + ^J n - a {kR)J n+1 _ a {kR) = Q for n = l,...,oo, (31b) 

when a > 1/2 (k > 0) while, for a < 1/2 (k < 0), the first equation in (|3l| ) holds for n — 1, ...,oo and the second one 
applies for n = 0, ...,oo. 

Of course, a regularization method must be introduced in order to give sense to the divergent sum in Eq. (|30|). In 
the framework of the ^-regularization PTlpq] (for several applications see P9LpO[ ) , 

Ec = ~M lim M 2 X(s) = ~M lim i 2^ (^^-) , (32) 

where the parameter M was introduced for dimensional reasons. 



Here, it is useful to define the so-called partial zeta function 

oo 

where k^j are the roots of 

2to 



J'(kR) - J'(kR) - —J^kRVn-iikR) = 0. 



(33) 



(34) 



So, after introducing v = n+ \ and a = a — i, the Casimir energy for the behaviour I at the origin can be written, 
for any a, as 



Eh = ~M lim M 2 * J £ [C, +a ( S ) + C„-a(*)] " Cj-|a|(«) 1 . 

U=i,i,... " J 



(35) 



while for the behaviour II at the origin, it is given by 



E, 



ii 
c 



l -M h mi M 2S J ^ [C +Q ( S ) + C„-a(*)] - Cl-sgn(«)a(s) 1 

U=i,|,.. J 



(36) 



From Eq. (Bq) it is clear that, as mentioned before, the Casimir energy for a behaviour of type I at the origin 
is independent of the integer part i of the reduced flux. Moreover, it is invariant under a — > —a (a — > 1 — a). 
Thus, it is enough to study it for < a < g, where the absolute value in the last term can be ignored, and to use 

£&(<*) = £&(-<*). 

Similarly, from Eq. (|36|), the Casimir energy for behaviour II is seen to be invariant under k — > — k. Thus we only 
study the case n > 0, where the last term is again (i_ a (s) and — i < a < ^ is considered. 



V. EVALUATION OF THE CASIMIR ENERGY 



The Casimir energies in Eqs. ( |35| ) and (R6J) contain two contributions: the term inside the square brackets, which 
is summed over v, and the last term, which is a partial zeta function. 

In both cases, it is useful to introduce, as in references |26]fil| , |4l| , an integral representation for the partial zeta 
function: 



Us) = 2 R 2i 



dx [x 2 -z 2 ]~ s ^- In x-^-^F^x) 



dx 



where 



F^x) = I 2 Jx) + I 2 x {x) + ^I^I^xix) , 



(37) 



(38) 



which have to be summed according to Eqs. (|35]) and (|36|). Here, the dimensionless variable z = mR has been 
introduced. 

In order to identify the divergences and evaluate the finite parts of the terms in Eqs. (|3g) and ( |3q ) an analytical 
continuation of the zeta function £(s) to s = —\ has to be constructed. A method to do this has been developed 
in J43J and for details of the procedure we refer to this reference. For the part of the zeta functions involving the 
angular momentum sum the method consists of adding and subtracting several orders of the uniform Debye expansion 
of Eq. ([58|) such as to make the sum as well as the integral in Eqs. (|3^), (36) and (|37]) well defined in an increasing 
strip of the complex s-plane. For the partial zeta function, instead, subtracting and adding the asymptotic terms for 
large arguments of the Bessel functions will be sufficient. 

Let us first study the terms summed over v. By making use of the recurrence relations for Bessel functions, it is 
immediate to obtain 



T(n,x,z) = —In 



x -2(p-l) F ^ x) 



2 (jJL 



(-)- 



+ *- (£) 2 * + 2^*) + £(£) 2 *<$(*) 



^ W 1 + (f ) 2 + I (f ) 2 - + I (f ) 2 dp (*) + ^ (f ) 2 ^,W + 7^ (f ) 2 <?W ' 



(39) 



where d fl (x) = x-p ln/ M (x). This expression can be developed in powers of —, through the Debye expansion of Bessel 



functions, after taking (^j 



Vl-7 



with r the variable of the recursive polynomials Uk (r) |45| 



If D^ N '(/j,, x,t, z) is such an expansion up to the order -^-, the partial zeta function can be written as 



Cm 00 = C(s) + tf(*) : 



where 



COO 



-R z 



dx \x 2 — z 2 ~\ 



T(n,x, z) — Z)' '([i, x, r, z) 



is the analytic part of the partial zeta function for s = — A , while 



Cfa) = 2 



SU17TS 



-R 



2« 



dx [x 2 -z 2 ] s £> (Ar V,x,r,z) 



(40) 



(41) 



(42) 



is the asymptotic contribution. 

In order to make the integral in Eq. (|4l|), and the subsequent sum over v : absolutely convergent at s = — A, it is 
necessary to take N > 2 p3| . We will choose N — 4 to improve the convergence of the sum of the analytic term (Ell) , 
thus decreasing the computational time needed to get accurate numerical results. 

Now, the term in square brackets in Eqs. (35) and ( pq ) involves the combination C, v+OL [s) + C»-a(s)- In order to use 

expansions, 



we further expand in powers of -. We introduce the corresponding combination of asymptotic 



2at 2 



v 
2at 2 



AW(i/ )I ,t,«) = DWUl + 2) lI| t(l + 2) 

expanded up to the order -i-. In the above expression £ = -= — = 
The asymptotic expansion can be written as 

N 

A( N V, x, M) = A-i + Ao + ]T A 2 , 
where 



( 1+ 2^ 



,z + 



27 



(43) 



i=i 



A_i 



4x t 

7l + t' 



An = 



2x £ 2 
^2 1 + i' 



1 rf v-^ 

A 4 = — — 2^ ft (i,i)* 



'+7 



^' dx 



(44) 



(45) 



i=o 



and the coefficients buj) are listed in Appendix |A|. 

Then, the term inside the square brackets in Eqs. (p3) and (Bq) can be written as 



oo 1 / N-l 

Urn M 2s £ [C+a(s) + C*-a(*)] = ^(-g) + *™ Af 2s A_!( S ) + A,(s) + £ ^(s) 



where 



Z{ -2 



irMR 



E 

ra=0 



r/.r 



[x 2 - (z) 2 ] 2 JT (v + a, x, z) + T (i/ - a, x, z) - A (JV) (i/, x, i, z)} 



(46) 



(47) 



A _ l{s) = 2 ^i R ^j2,rd X [(, x r- z r s ^^^ 1 - 

n=0 ~ 



(48) 



A ( s ) = 2 S -^^R 2s 2 V / dx [{vx) 2 - z 2 ] 



n=a 



1 Vl + a; 2 - 1 



VTT^ 2 



(49) 



Ai (s) = 2^R 2 °'£b m J2v- i I dx [{vxf - z 2 ]~ s ^ 



1 



rfx \y/l 



, p z 



<+J 



j=0 »=0 

Eqs. (f48|), ( ff9| ) and ( J50J ) can be expressed in a systematic way by introducing the functions |25| 

, 2" 



/(*;o l 6;a;) = £i^ 1 + g) 



n=0 



-s-b 



studied further in Appendix pi which allow one to write the asymptotic parts as 



(50) 



(51) 



(52) 



*M-3W— iX^'M' 



(53) 



2/ 



where 



3=0 



^ (i)i) ( S ) = - 2 ^^-(^) Ife|l ^ 2 7 f s;i , i±i. ti 



r(^)r( s ) 



(54) 



(55) 



The complete expressions for these asymptotic parts around s = — i are derived in Appendix |C[ Here, we list their 
residues, which will be relevant to the discussion of the renormalization in the next section: 



ResL _i A-x =0, 



Res| s= _i A = - 
Res\ s= _iAx = 0, 



ResL_ i Ai = 



R 






1 1 a 2 z z z 2 

64 12tt ~~ ~ 4 7T Y 



(56a) 
(56b) 
(56c) 

(56d) 



Next, we study the partial-zeta contribution 



1 



e c = --M lim M 
2 s ^- 



2 s 



"&-«(*) 



(57) 



to the Casimir energy in Eqs. ( p5| ) and (|3^), and for reasons explained we omit the absolute value in the index. 

In order to isolate the singularities, it is enough to consider the three leading terms in the asymptotic expansion of 
Bessel functions for large arguments, which will be denoted by L(-| — a, x, z) ; thus, the partial zeta function can be 
written as 



Ci-«(*) = <f_ a (s)+d a ( s ), 



(58) 



where 



C? _ ( s) = 2 — i? 2s / dx \x 2 -z 2 Y S 

2 7T Jz 



T I 2 ~ a ' x ' z J ~ L ( 2 ~ a ' X ' Z 



(59) 



C rf_ a ( s) = 2 ^li? 2s /""(fa [^-^""Lfi-a.a:,* 



(60) 



Now, the subtracted terms can be written as 



L ( - — a, x,z ] = <5q + #i + fo • 



(61) 



where 



<5 = 2, 5i = -2a, <5 2 = -^{a 2 - z) . 



(62) 



Then 



lim M 2s (i_ a {s) = z(--)+ lira M 2s (a (s) + oi(«) + o 2 («)) 

>^— i 2 2 s— > — i 



(63) 



with 



{ - 1 2) = -Jir[ ^[x"-z 2 ]nT(\-a,x,z)-L(\-a,x,z 



(64) 



which will be evaluated numerically, and 



a (s) 



R 



z 2 1 



1 + 21 ° g (^) 



O s 



(65) 



Ol(«) = R 



2az + O [ s + - 



(66) 



a 2 (s) = —(a 2 - z) { - p- + - 



2 + 21og(-M 



\2RJ 



+ 0[s + 



(67) 



So, the residues at s = — o are given by 



Res| s= _ia = - 
Res| s= _ i ai = 0, 

Res| s= _i a 2 = - 



(68a) 
(68b) 

(68c) 



VI. DISCUSSION OF THE RESULTS 



Clearly the Casimir energy is divergent and using Eqs. ( |56| ) and Eqs. (Bq) in Eqs. (pq) and 
is given by 



the total residue 



Res| s= i£ c 



1(1 z z 1 



2R 64 4 2 



(69) 
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FIG. 1. Difference Ed of Casimir energies. Behaviour I at the origin. Top to bottom: z = yig, |, |, |, i. 



which is independent of the flux. Thus, the difference between Casimir energies with arbitrary and integer flux is 
finite and contains the relevant information about the effect of the flux. 

In Fig. [j], we plot the dimensionless difference Ed = R (Ec(a) — Ec(0)) for a behaviour of type I at the origin, as a 
function of a (the fractionary part of the reduced flux), for different values of z. Since the finite part of the Casimir 
energy is continuous in a, the difference goes to zero both at a — and a = 1. It shows a minimum at a = -. It 
is interesting to note, that around a = 0, the vacuum energy decreases when the flux grows and that this effect is 
enhanced with increasing mass. 

For a type II behaviour at the origin, the same difference is plotted in Fig. @. With decreasing mass our curves tend 
to the m — result of Ref. |2q| and already z — -r^ shows quite a good agreement with the corresponding figure in 
that reference (except for a factor of 2, due to the fact that only one polarization was considered by the authors of 

{§)■ 

Whereas for small values of the mass the energy exhibits a minimum at a ^ 0, for larger values of m this minimum 
is shifted towards a = 0. Furthermore, for m / 0a nonzero value is seen to arise for a — > + . This is due to the 
discontinuous behaviour of the finite part of the vacuum energy, more precisely of the contribution of the critical 
subspace, at integer values of the flux. 

The origin of this discontinuity can be traced back to the appearance, for a — ► + , of a root of the combination of 
Bessel functions involved in the partial zeta ( a - Such a root is absent when a = 0. For a — ► + , this root goes to zero 
and, thus, gives rise to a gap, which equals m. The quantity J^(kR) — J a l _ 1 (kR) — ^J a (kR)J a -i(kR) is shown in 
Fig. H as a function of kR, for various values of a in order to clarify this discontinuous behaviour. 

In summary, we have seen that the presence of the mass as well as the choice of the self-adjoint extension of the 
Hamiltonian have a considerable influence on the dependence of the Casimir energy of the flux. We have analyzed and 
discussed in detail the behavior of the Casimir energy as a function of the parameters, namely flux and mass. As a 
consequence of the mass, in the case of the self-adjoint extension, (13), a discontinuity of the energy is found. It would 
be interesting to see how this effect comes about, starting with the finite radius flux p] and shrinking the radius to 
zero. Does the effect persist or it is only a result of the singular vortex? Furthermore, one should try to understand 
better the physical meaning of the different self-adjoint extensions by considering how the Casimir energy depends 
on the parameter of the one-parameter family of self-adjoint extensions. Finally, more realistic (3 + l)-dimensional 
calculation should be envisaged. 
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APPENDIX A: COEFFICIENTS 



In this Appendix we list the coefficients buj-j defined by Eq. ( (45| ) 
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APPENDIX B: FUNCTIONS f(s;a,b;x) 

Here we are going to provide all analytical properties of the functions f(s;a,b;x) defined in Eq. (pl|). 
As in Rcf. pq |, we will make use of 



„=o J o jo l + e 



in the limit e — > 0. 

When applied to h{v) = v a 1 1 + (-) j , the previous equation gives 



n=0 



na\ j u r ^ 2\~ t 



r 1 u a 
2sin( — — ) / du 7, (l — u 

POO 

/7ra \ I u 



2sin^-^ j / i du—^-^-ipj. (B2) 

Now, we are interested in f(s;a,b;x), for arbitrary b. F rom the definition, it is clear that f(s;a,b;x) = f(s 



b,a,0;x). However, as b grows, the integrals in Eq. (B2) eventually diverge at u = 1. In order to avoid such 



divergences, we will perform an adequate number of integrations by parts, thus obtaining 
f(t;a,0;x) = 



/+ Jirmr(t-^; 



2 T(t) 

e" a + l\ (-l)f r ( f _a) /-oo 

j^-sm(Trt) V 2 > j udug(-,a-l;u,x)(u 2 -l)- t+ - 



2 

.i\ 1 rf*.tlUfi / a _i \ 

udug I — - — , a — 1; w, x I (1 — u ) 



e i7ra -l\ 1 r(t-2=i) r Z" 1 . /" I 





a- 1 



2^ r(t) 

(— l) 2 * - cos(7ri) / udug ( , a — 1; u, x I (w 2 — l)~ t+ ~ 



where 



(B3) 



g(a,b;u,x) = [ -j-j-) : , " . (B4) 



However, the number of integrations by parts is bounded by the requirement that the integrated terms are well- 
behaved at u = 0. In what follows, we will thus keep the number of integrations admissible, by making use of the 
following recurrence relationship 

f(s; a, b; x) = f(s; a,b-l;x) j/(s; a + 2, b; x) . (B5) 

x 

In this way, all the required functions can be reduced to four different cases 

• f(s; 2n, n; x), f(s; 2n,n + h;x), n = 0, 1, 2, 3, 4, 5, 6, 7 

• f(s; 2n + 1, n; x), f(s; 2n + 1, n+ ~; x), n = 0, 1, 2, 3, 4, 5, 6 . 
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Finally, after expanding in powers of s + i, we get 
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(B7) 
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[1 + n + n [-0(1) - ip(n + 1)]] + 



(B9) 



where ip(x) is the Euler Psi function. 

These expressions generate all the /-functions necessary for the evaluation of the required Ai (s), for i = 1, 2, .... 

Notice that, for j = and i = 1, the prefactor T (s + ^-) in Eq. ( p5[ ) has a pole at s = —\. Thus, the order s + 5 
in the expansion of /(s; 0, | ; x) must be retained, 



/ I s;0, -;x 



-7TX S + 



1 +2 / udug(Q,—l;u,x) 



'rOiiS+tf 



(BIO) 



APPENDIX C: EVALUATION OF A-i AND A 

Finally, in this Appendix we are going to describe some details of the calculation of the Ai(s), Eqs. (p2|)-(p5|), at 



We will, in the first place, obtain an expression for A_i in equation ( |52| ) 
By using equation (|B2J ) in the previous appendix, it can be put in the form 

u ; V^ r(«) i Vy \2 V ' ,/ r( s -|) 



(Cl) 



2v^zsin tt(- - s) / dw (u 2 - l) 



1 



1 + e 2 ™^ 2 



After interchanging the integrals, one gets 



(C2) 
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A-i(s) 
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oo x n 



where Li,- (a) = E„=i |r 

Finally, expanding around s = — i, one has 
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which is a useful representation for numerical calculations. 
Let us now go to the evaluation of Aq in Eq. (J53|) , 
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/? r(s) 
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As before, using Eq. (B3), and interchanging the integration order, one gets 
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Now, after analytically extending, and developing around s = — -|, one has 
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7Ti?( S +i) V 12 
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27TU2 
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where a simple pole at s = — h appears. 

Clearly, in both cases the finite parts must be evaluated numerically. We will not go into the detailed calculation 
of the Ai for i > 0, since it is a direct consequence of the properties of /(s; a, b; x) described in the previous appendix. 
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